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The related studies of two-particle correlations and event-by-event fluctuations have played im- 
portant roles in the search for new physics through the experimental study of relativistic heavy ion 
collisions. We present a general method of determining minimum-bias two-particle correlations 
and show the relationship between these correlations and event-by-event fluctuations. Data from 
the STAR experiment at RHIC for Au-Au collisions at a/Saw = 62 and 200 GeV will be presented 
that show the energy and centrality dependences of angular correlations. This analysis provides a 
unique method for studying the interaction of semi-hard scattered partons with the dense medium 
produced at RHIC as well as properties of the bulk medium itself, and will shed new light on the 
sources of non-statistical fluctuations. 
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1. Introduction 

The analysis of correlations and fluctuations provides essential information on the nature of 
the medium produced in ultra-relativistic heavy-ion collisions. For example, correlation studies 
often focus on specific effects such as the collective behavior leading to elliptic flow, quantum 
correlations, and high-/?, jet correlations. Event-by-event fluctuation analyses of global parameters 
(e.g. mean p t , net charge, etc.) search for critical fluctuations as a signature of a phase transition. 
This dichotomy is evidence of a false perception that correlations and fluctuations are unrelated. 
This work presents a general method of measuring minimum-bias two-particle correlations, where 
a two-particle correlation is constructed using all detected particles from an event combined into all 
possible unique pairs, that includes the above mentioned individualized correlation analyses and is 
directly related to event-by-event fluctuations. 

In Section 2 a correlation measure is derived from standard statistical definitions. In Section 
3 this method is used to measure angular (pseudorapidity r\ and azimuth 0) correlations in 62 and 
200 GeV Au-Au collision data from the STAR detector at RHIC [JXJ] _ These results, combined 
with similar analyses at 130 GeV and of 200 GeV proton-proton collisions, present a survey of 
the system energy and centrality dependence of two-particle correlations at RHIC. The evolution 
with centrality and comparison to proton-proton events and HIJING simulations indicate that par- 
ton scattering and fragmentation are modified in central heavy-ion collisions, presumably through 
interactions with a dense medium. In Section 4 the technique is extended to differentiate between 
electric charge sign to form charge-dependent correlations which provide differential information 
on hadronization in contrast to that inferred from net-charge fluctuations. 

2. Correlation measures 

Two-particle correlation measures are formed by comparing the two-particle density p(p\,P2) 
to an appropriate reference scaled by a normalization factor. The object of interest is the two- 
particle density of pairs of particles chosen from the same event, referred to as sibling pairs, defined 
as PsibiPiiPi)- Pairs of particles from different but similar events provide a mixed-event reference 
Pref{P\iPi) containing important detector acceptance and phase space effects. These two factors 
may be combined in a number of ways. For analysis of quantum correlations (HBT) p s $ / p re f 
is projected onto a small region of relative momentum space where approximately every pair is 
correlated. Another common choice is Ps,b p ^ nf , where the ratio removes detector acceptance and 
provides cancellation of systematic effects. These ratios are per pair measures and thus contain 
a 1/JV dilution factor for average multiplicity N. Measures containing such dilution factors are 
not directly comparable across system energies, centralities, incident particle species, or any other 
condition affecting multiplicity, but more importantly per pair measures can be misleading. For 
example, elliptic flow measure V2 includes a N~ 1 / 2 scaling, giving the impression that elliptic flow 
is quite large in peripheral collisions, and therefore does not smoothly extrapolate to the proton- 
proton limit where the cos[2(0i — fc)] correlation component is negligible. 

Two-particle correlations are often measured in the context of a high-/?, triggered and associ- 
ated particle, leading to a model-dependent jet hypothesis and trigger bias, however this condition 
is not necessary. The same analysis can be performed constructing every possible pair for each 
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event (i.e. "triggering" on every particle), creating a minimum-bias correlation having several ad- 
vantages over a conditional correlation. As particles interact with the medium in RHIC collisions 
they lose energy, placing particles which have significant interactions potentially below the trig- 
gered momentum region. Removing this restriction allows the search for particles with enough 
momentum to be easily detected, yet soft enough to have significant interactions with the medium, 
conceptually similar to finding an ideal Brownian probe for heavy-ion collisions j2|]. It is also pos- 
sible to use correlations to explore the bulk medium itself, such as the response of the medium to 
a hard scattered parton pp.. Finally, most correlation analyses must extract the correlation signal 
of interest, such as jets, HBT, or flow, from other processes, but a minimum-bias analysis gives a 
complete picture of the dynamics of the system and far more information to use in determining the 
relative contributions from these components. 

It is possible to derive a correlation measure from standard statistical definitions that retains 
the useful properties of removing detector acceptance effects and the cancellation of tracking in- 
efficiencies and other systematic errors in the ratio p s ib/Pref, but also scales with N in way that 
allows for direct comparison between different collision systems, energies, and centralities. One 
way to achieve this scaling is to create a correlation measure that would be independent of central- 
ity if heavy-ion collisions were linear superpositions of proton-proton collisions, making deviations 
from this participant number scaling directly accessible. The canonical definition of a correlation 
between two variables, usually referred to as Pearson's normalized covariance or correlation coef- 
ficient, is the ratio of the covariance between the two variables with the product of their standard 
deviations: corr(a,b) = co ^^ ■ This work will focus on correlation of the number of pairs of 
particles, but may easily be extended to other quantities. Event-by-event fluctuation measures can 
be defined in terms of an integral of corr(a,b), but examining the integrand directly reveals new 
differential information of the system. Only correlations will be shown here, but the relationship 
between fluctuations and correlations is emphasized to place these results in a larger context and 
highlight the applicability of this analysis method. A related study [||] illustrates this by explicitly 
inverting a transverse momentum fluctuation integral equation to find the associated correlations. 

Returning to the definition corr(a,b) = Efngp f we must express densities p sl j and p re f as a 
covariance. Particles from different events must be uncorrected, so the two-particle density p re / 
can be factorized as a product of one-particle densities. Defining Ap = p^ — p re f we see that in 
terms of particle counts n in histogram bin (a,b) (e a e^)Ap = n a rib — n a nt, = (n — h) a (n — n)b, where 
£ is the bin width used to convert densities to bin counts and the overline denotes an average over 
events. Thus Ap is a covariance, which is more conveniently expressed as particle counts than 
densities. In the Poisson limit o a Ob = ^n a n~b, defining our correlation measure as the density ratio: 



Ap = 1 [n-n) a {n-n) b 



£ a £b 



yjn a n h 



The correlation measures Ap, Ap/p re f, and Ap / yjp re f have very similar structures in an- 
gular (pseudorapidity v\ and azimuth <p) space when corrected for acceptance so that p re f is ap- 
proximately uniform, though with different overall normalizations. Of these, only Ap / ^/p re f is a 
per-particle measure which allows for direct comparison of correlations from events regardless of 
multiplicity. The square root in the denominator provides a bin-by-bin normalization rather than 
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applying a global normalization of some factor of N across all bins. This distinction is vital in 
transverse momentum space where both p s $ and p re f are rapidly varying. 

3. Experiment 

Data for this analysis were obtained from the STAR detector at RHIC measuring charged 
particle production from Au-Au collisions at ^snn of 62 and 200 GeV. A minimum-bias event 
sample required coincidence of two Zero-Degree Calorimeters, a minimum threshold in the Cen- 
tral Trigger Barrel, and a reconstructed primary vertex. Particle tracks were reconstructed with the 
Time Projection Chamber operating within a uniform 0.5 T axial magnetic field. Tracks satisfying 
standard quality requirements were accepted within full azimuth, |t}| < 1.0, and 0.15 < p t < 6.0 
GeV/c. Electrons and positrons were suppressed with a dE/dx cut within 0.2 < p < 0.4 GeV/c, 
but no other particle identification was attempted. Pair-wise tracking inefficiencies were corrected 
by applying track merging and splitting cuts to both sibling and mixed-track pairs. To reduce sys- 
tematic error, events were grouped into multiplicity and primary z-vertex classes such that within 
a class AN < 50 and Az, < 5 cm. Reference pairs were created by only mixing events within the 
same class. Furthermore, the correlation measure was constructed for each electric charge pair 

type (i.e. ++, H — , — h, and ) for each class, and then combined within each centrality bin. 

Events for both energies were divided into 1 1 centrality bins defined by total cross section fraction 
as 90-100% (peripheral), 80-90%, 70-80%, . . . , 10-20%, 5-10%, and 0-5% (central). 

The two-particle correlation was measured as a function of angular difference variables 
77a = T]i — T]2 and A = 0i — fa. Correlation structures are found to be invariant along the sum 
variables Tji + rj 2 and 0! + 2 , thus projections onto coordinates (t]a, 0a) compactly represent the 
entire 4D angular space (r\\ , r/2, 0i , 02) without loss of information. 

4. Charge-independent correlations 

To understand two-particle correlations in the complex system of heavy ion collisions, we 
first examine proton-proton (p-p) collisions as a reference, p-p angular correlation data are well 
described by a two-component model of string plus jet fragmentation, a similar model is often 
used by event generators such as PYTHIA [§], HIJING and AMPT [§] to simulate particle 
production mechanisms. Analysis of the multiplicity dependence of transverse momentum spectra 
in p-p collisions showed that the spectra can be separated into soft and semi-hard components [f7|]. 
Based on this result, a follow-up study of two-particle correlations observed that the soft and semi- 
hard components were distinctly separated in transverse rapidity y t 1 where the soft component was 
localized at y t < 2 (p t < 0.5 GeV) for each particle in the pair, while a Gaussian distribution on y t 
peaked at y t ~ 2.7 describes the semi-hard region 

Using y t as a cut space, the angular correlations for soft and semi-hard components were ex- 
amined. The semi-hard component correlations produce a large peak centered at the (tja, 0a) origin 
plus an away-side ridge at 0a = % and are fully consistent with conventional high-p f jet angular 

'Transverse rapidity y t = | ln( ) = ln( m '^' ) with transverse momentum p,, transverse mass m, = ' pj + itIq 
and mass mo (pion mass is assumed for unidentified particles). It is quite natural to describe transverse fragmentation 
with respect to transverse rapidity in analogy to analysis of longitudinal string fragmentation. 



4 



Angular Correlations in STAR 



Michael Daugherity for the STAR Collaboration 



correlations, even though the measurements extend to a much lower momentum range and use no 
leading particle trigger. The observation of jet-like correlations at a lower momentum range under- 
pins the terminology of semi-haid scattering and minijet (minimum-bias jet) to distinguish these 
correlations from those obtained by conventional high-/?, triggered particle analysis. One surpris- 
ing result is that minijet correlations in p-p collisions are not symmetric on relative pseudorapidity 
and azimuth. Instead, the width on 0a is nearly 50% larger than the width on tja, in strong disagree- 
ment with an angularly symmetric jet cone. The soft component contains a large HBT peak (only 
observed in like-sign pairs) at the origin and a ID Gaussian centered at 7]a = that is independent 
of 0a> consistent with longitudinal string fragmentation. This two-component model provides a 
framework for understanding correlations in RHIC heavy-ion collisions where the same partonic 
scattering processes observed in p-p are assumed to be embedded in an unknown medium. 

A previous analysis reported angular correlations in Au-Au collisions at 130 GeV (us- 
ing correlation measure N-r^-) showing substantial differences between p-p and Au-Au as well as 

Pref 

substantial evolution of observed correlation structures with centrality. That analysis accepted a 
transverse momentum range of 0. 15 < p t < 2.0 GeV/c, containing both soft and semi-hard compo- 
nents. The soft component, represented by a ID Gaussian at tja = 0, is diminished relative to p-p 
in the peripheral centrality bin and is not observed in mid-centrality, suggesting that longitudinal 
string fragmentation is not a relevant mode of particle production in central heavy-ion collisions. 
The semi-hard component contains a minijet cone that in p-p is elongated on 0a » and in Au-Au 
deforms with increasing centrality to become significantly more broad on T]a and slightly narrower 
on 0a> reversing the asymmetry found in p-p. The most straightforward interpretation of the elon- 
gation in tja is to assume that the same semi-hard partonic scattering processes measured in p-p 
collisions are now embedded in the longitudinally-expanding medium created in heavy-ion colli- 
sions. Multiple hadronic scattering and other final-state interactions would most likely broaden the 
minijet peak along both tja and (j) A , and are unlikely to be of sufficient magnitude to explain this 
effect. However, partonic interactions with a longitudinally-expanding medium could be the source 
of the observed asymmetry reversal. 

Unfortunately, the limited volume of the 130 GeV Au-Au data only allowed four centrality 
bins, the most peripheral being 40-70% of the total cross section, and could not reveal the full 
centrality dependence of correlation structures and the transition in peripheral collisions from p-p 
to Au-Au. The present study analyzes nearly an order of magnitude more events at both 62 and 200 
GeV which are divided into eleven centrality bins spanning the entire range from 0-100% of total 
cross section, revealing for the first time the detailed centrality and energy dependences of these 
correlation structures, as shown in Figures |and[| The data span the transverse momentum range 
0.15 < p t < 6.0 GeV/c and, as in the 130 GeV analysis, contain no trigger particle of correlation 
model assumption. 

These correlation structures may be summarized by three general trends. First, the most pe- 
ripheral (90-100%) centrality bins are in excellent agreement with p-p data, and the correlation 
structures evolve smoothly with increasing centrality for each energy. Second, the soft component 
(ID Gaussian at tja = 0) diminishes with increasing centrality and does so much more rapidly 
with increasing energy. The soft component is visible at 62 GeV until ~20% centrality whereas 
at 200 GeV it is completely absent by 70%. Finally, the hard component minijet peak is dramati- 
cally broadened in tja with increasing centrality and this deformation occurs much more rapidly at 
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higher energy. 

The HIJING Monte Carlo simulation model includes particle production from soft fragmenting 
longitudinal strings and hard perturbative QCD jet production and fragmentation with jet quench- 
ing. The results from HIJING simulations qualitatively agree with the peripheral data, though 
HIJING substantially over-predicts the amount of soft component relative to hard and fails to pre- 
dict the energy and centrality dependence. HIJING also fails to predict the asymmetry of and 
0a 5 which is unaffected by the jet quenching option, and actually predicts very little change with 
centrality. These results suggest that the observed broadening along may be due to microscopic 
interactions with the medium not modeled by HIJING. Future studies including partonic interac- 
tions within the AMPT model will examine this effect in more detail. 



5. Charge-dependent correlations 

The above analysis can be extended to differentiate between the relative electric charge sign 

of particles. Defining like-sign pairs (LS: ++, ) and unlike-sign pairs (US: H — , — h) the above 

charge-independent (CI) correlation is formed as CI = LS + US. A charge-dependent (CD) cor- 
relation is defined as CD = LS - US, where the sign convention is chosen to be compatible with 
isospin conventions 2 . Charge-dependent angular correlations in p-p collisions are dominated by a 



negative ID Gaussian on tja with width oVja ~ 1 [10, 11]. This feature is conventionally attributed 



to charge-ordering during longitudinal string fragmentation Q12| ] and shows the dependence of the 
soft component on relative charge sign. Since HBT correlations are only present in LS pairs, they 
form the only other significant correlation feature of a positive Gaussian peak centered at the origin. 

In central Au-Au collisions at 130 GeV the ID Gaussian is not detected, instead a negative 2D 
exponential dominates the correlation structure, suggesting charge-ordering is no longer occurring 



along a ID string as in p-p but along a 2D surface [ 13Q. The change from Gaussian to exponential 
distributions may be caused by pair attenuation which increases with opening angle due to inter- 
actions with a the medium. Lower centralities show intermediate stages between p-p and central 
Au-Au, but again this analysis only shows coarse centrality dependence. 

The present analysis shows angular CD correlations in Au-Au collisions at 62 and 200 GeV in 
Figures || and |j. As in the previous section, excellent agreement between the most peripheral (90- 
100%) centrality bin and p-p is demonstrated and the correlation structure evolves smoothly with 
increasing centrality. These results confirm the independent observation at 130 GeV that the ID 
Gaussian on tja in p-p, attributed to longitudinal charge-ordering, is replaced by a 2D exponential 
approximately symmetric on T7a and A - The evolution occurs more rapidly at higher energy with 
a more dramatic decrease in widths of the exponential. 

The primary source of CD correlations appears to be charge-ordering due to local charge con- 
servation in fragmentation processes into quark-antiquark pairs. It is more likely for neighboring 
particles to have to opposite charge, producing a stronger correlation in US than LS, thus the CD 
correlation is negative. The symmetry axis of charge-ordering is telling of the fragmentation pro- 
cess. In elementary collisions charge-ordering can be seen on the longitudinal (z or rj) direction 



2 The isovector term in the two-particle density is the coefficient of Ti • T2 where T is the z-component of isospin. 
For a pair of LS charged pions X\ ■ T2 = +1 while US Tj • T2 = — 1, thus CI = LS + US isolates the isoscalar component, 
while CD = LS - US is isovector. 
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from string fragmentation phenomenologically described by the Lund model [12]. It is also possi- 
ble to observe transverse charge-ordering associated with jet fragmentation processes on transverse 
rapidity, which are strikingly similar to CD correlations of soft longitudinal fragmentation on lon- 
gitudinal rapidity. 

In Au-Au collisions at RHIC, we observe a significant evolution of the charge-ordering signal. 
Peripheral CD correlations, which are in excellent agreement with p-p data, show a ID charge- 
ordering signal on tja as expected from fragmentation of ID longitudinal strings. HIJING predicts 
no change in CD correlations with energy or centrality, but the data show otherwise. We observe 
that ID charge-ordering on rj^ evolves to become 2D charge-ordering on tja and 0a> measuring/or 
the first time charge-ordering in the azimuthal direction. These results suggest that in central Au- 
Au collisions particle production does not occur along several independent ID longitudinal strings, 
but along a surface in at least two dimensions (the transverse direction is not considered here). 
Additionally, we also find a change in shape with increased centrality of the charge-ordering signal 
from Gaussian to exponential. This change is expected in a scenario where rescattering increases 
with the opening angle of the pair. While hadronic rescattering may contribute to this effect, it 
is extremely unlikely that rescattering alone could cause a ID signal to appear to be both 2D and 
symmetric without invoking new hadronization mechanisms. 

Several net-charge fluctuation measures have been proposed and studied to search for QGP 
formation [14]. Those analyses have shown little or no centrality dependence in sharp contrast 
to the predictions of dramatic suppression of net-charge fluctuations. We must again emphasize 
that net-charge fluctuation measures typically integrate over the correlations shown here and are 
therefore insensitive to changes in differential structure. It is significant to note that while we 
observe a dramatic change in the amplitude, width, and shape of charge-ordering correlations with 
centrality, the integrals of these correlations exhibit only modest changes with centrality (at 130 
GeV the integrals increase by about 20% from peripheral to central [13]). 



6. Summary 

The primary goal of the heavy-ion physics program at RHIC is to determine if the nuclear 
matter created in RHIC collisions is of sufficient energy density such that the relevant degrees of 
freedom are partonic rather than hadronic and, if this partonic medium is created, to measure its 
properties. To that end we have presented a general analysis technique based on minimum-bias 
two-particle correlations which are shown to be sensitive to the low Q 2 kinematic range where 
scattered partons strongly interact with the bulk medium. In addition to providing access to new 
physics, this analysis represents a generalization of both triggered jet correlations and event-by- 
event fluctuation analyses, enabling results from these methods to be evaluated in a much larger 
context and from the point-of-view of a more complete physical picture. 

Surveying the correlations in Au-Au collisions at 62 and 200 GeV, along with previous anal- 
yses at 130 GeV and in proton-proton collisions, provides an overview of the centrality and en- 
ergy dependence of correlation structures from heavy-ion to nucleon-nucleon collisions. Charge- 
independent correlations reveal the attenuation of soft longitudinal string fragments and dramatic 
broadening of the minijet correlation peak, both phenomena are consistent with interactions with 
a partonic medium developing with increasing centrality and collision energy. Charge-dependent 
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correlations suggest a fundamental change in hadronization geometry from one dimensional longi- 
tudinal strings to a multi-dimensional surface, also consistent with the hadronization of a partonic 
medium. 

I thank the organizers of the 2nd International Workshop on Correlations and Fluctuations 
in Relativistic Nuclear Collisions. This research in supported in part by the U.S. Department of 
Energy. 
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Figure 3: Au-Au 62 GeV CD angular correlations. Note changes in scale. See text for discussion. 
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Figure 4: Au-Au 200 GeV CD correlations. 



10 



